Mechanical effects of optical resonators on driven trapped atoms: 
Ground state cooling in a high finesse cavity 
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We investigate theoretically the mechanical effects of light on atoms trapped by an external 
potential, whose dipole transition couples to the mode of an optical resonator and is driven by a 
laser. We derive an analytical expression for the quantum center-of-mass dynamics, which is valid in 
presence of a tight external potential. This equation has broad validity and allows for a transparent 
interpretation of the individual scattering processes leading to cooling. We show that the dynamics 
are a competition of the mechanical effects of the cavity and of the laser photons, which may 
mutually interfere. We focus onto the good-cavity limit and identify novel cooling schemes, which 
are based on quantum interference effects and lead to efficient ground state cooling in experimentally 
accessible parameter regimes. 



I. INTRODUCTION 

Atom cooling by photon scattering is achieved by en- 
hancing the rate of scattering processes that dissipate 
motional energy, thereby exploiting the conservation of 
internal and mechanical energy in the interaction be- 
tween atoms and electromagnetic field [ljj. The atomic 
scattering cross section can be significantly modified by 
the coupling to an optical resonator, which acts both on 
the internal as well as on the external degrees of freedom. 
Hence, the scattering properties can be tailored, allowing 
to achieve efficient cooling also for atoms and molecules 
which may not offer a convenient configuration in free 
space 0, H[ . This principle is at the basis of cooling by 
means of an optical resonator. Indeed, the mechanical 
effects on atoms coupled to an optical resonator are ob 
ject of several experimental 
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investigations, which aim at developing a systematic un- 
derstanding of these complex dynamics both for its fun- 
damental aspects, as well as for the perspective of a high 
degree of control of complex systems with scalable num- 
ber of degrees of freedom. 

In this work we investigate the cooling dynamics of 
atoms inside optical resonators, when their center-of- 
mass motion is tightly confined by an external po- 
tential, like for instance a dipole @, [H| or an ion 
trap [ll|, H3, HH, [26j]. We consider the situation where 
an atomic optical dipole transition is driven by a laser 
and by a cavity resonator, as sketched in Fig. [U and dis- 
cuss in detail the results presented in [2?| . In particular, 
we show the detailled derivation of the rate equation dis- 
cussed in [27j ■ This equation has broad validity, which is 
supported by numerical checks, and allows for a transpar- 
ent interpretation of the individual scattering processes 
leading to cooling. Moreover, in the corresponding pa- 
rameter regimes it reproduces the results reported in [131 ] 
and [3. 

In this manuscript we mostly focus onto the good- 
cavity limit. In this regime we discuss when efficient 



cooling into the potential ground state can be achieved. 
In particular, we show that in experimentally accessible 
parameter regimes one may obtain almost unit ground 
state occupation, even when the natural linewidth of the 
dipole transition would not allow for ground state cooling 
in free space. Efficient ground-state cooling is often found 
by exploiting interference effects, arising from phase cor- 
relation between the laser and the field scattered by the 
atom into the resonator. Most of these interference ef- 
fects are due to the discrete nature of the spectrum of the 
center-of-mass motion, which is trapped by a harmonic 
potential. Hence, the dynamics here studied differ sub- 
stantially from the ones of cooling of free atoms inside 
cavities [3 [It], 0, [H| • Such interference effects are 
at the basis of novel cooling schemes, some of which have 
been identified in [27[, and which are discussed in detail 
in the present work. 




FIG. 1: An atom is confined by an harmonic potential of 
frequency v inside an optical resonator. A mode of the res- 
onator couples with strength g to the dipole, which is driven 
transversally by a laser at Rabi frequency Q. The system dis- 
sipates by spontaneous emission of the atomic excited state 
at rate 7 and by cavity decay at rate k. The other parameters 
are discussed in Sec. Mil 

This article is organized as follows. In Sec. UT1 some pre- 
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liminary considerations are made. In Sec. IIIII the model 
is introduced, and the basic equations for the motion are 
obtained. In Sec. |lV]we discuss the dynamics of cooling 
from the rate equation we obtain and review previous re- 
sults presented in the literature. In Sec. fVl novel cooling 
schemes are presented, whose dynamics are due to quan- 
tum correlations which are established in the good cavity 
limit. In Sec. I VII the results are reported: The cooling ef- 
ficiencies in the various parameter regimes are discussed 
and compared. In Sec. I VIII the conclusions are drawn. 
The appendices report detailed calculations at the basis 
of the equations derived in Sees. IIIII IIVI andfVl 



dressed states 



II. 



MECHANICAL EFFECTS OF CAVITY AND 
LASER ON THE ATOMIC MOTION 
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FIG. 2: Sketch of the internal levels \g) and |e) of the atomic 
dipole transition, driven by a laser and a cavity mode with 
coupling strengths g and Q, respectively. The arrows show 
the cavity and laser frequency with respect to the dipole fre- 
quency. Here, A and 8 C are the detunings of atom and cavity, 
respectively, from the laser frequency, A c is the detuning of 
the atomic transition from the cavity frequency. The fre- 
quency of the atomic transition (loo), of the cavity mode (cj c ) 
and of the laser field (lol) are indicated in the vertical scale. 

In this section we make some physical considerations, 
in order to provide insight into the results presented in 
the rest of the manuscript. The scattering cross section 
of the bare atom is usually very informative about the 
cooling process [28| . When the atomic transition is driven 
in saturation, this analysis is more conveniently done in 
the dressed state picture [ltl [29|, [3(3]. At this purpose, 
we first assume that the atom is fixed at the position 
x, such that the coupling constant to the cavity mode 
is g = g(x). The dipole level scheme and the relevant 
parameters are shown in Fig. [21 We denote by \g, n c ) 
and e, n c ) the states of the system, where \g), |e) are the 
ground and excited state of the atomic dipole and \n c ) the 
number of photons of the cavity mode. In the situation 
in which the atom is strongly coupled to the cavity mode 
and weakly pumped by the laser, the states which are 
relevantly involved into the dynamics are |g,0 c ) and the 
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|+) = sin?9|g, l c ) + cos$|e,O c ) 
|— ) = cosz%, l c ) — sin^|e, C ) 



tan i? = g/(-A c /2 + V.9 2 + A*/4) 



(1) 

(2) 



and A c the detuning between cavity mode and atom. 
Setting the energy of \g, C ) at zero, the frequencies of 
the states |±) are 



A± = -A c /2± v /.g2 + A2/4 



(3) 



and the respective linewidths are 7 + ~ Ksin 2 $ + 7 cos 2 
7_ ~ « cos 2 1} + 7 sin 2 where 7 is the linewidth of the 
dipole transition and k the cavity decay rate. The weak 
laser probe couples the dressed states |g,0 c ) and |e,0 c ). 

Signatures of the dressed states are for instance the 
resonances in the rate of photon scattering as obtained 
by scanning the probe laser through atomic resonance. 
This situation is depicted in Fig. [3J Here, the curve has 
been evaluated for a good resonator, namely k < 7,3 and 
A c 7^ 0. For these parameters the linewidth of one of the 
two resonances is narrower than the natural linewidth of 
the dipole. Moreover, when the probe laser is resonant 
with the cavity mode, the spectrum exhibits a minimum, 
which reaches zero for k = 0, namely no photons are 
scattered. This behaviour is due to an interference effect 
between laser and cavity resonator, such that there is no 
radiation scattered by the atom, as it is at a point where 
the two fields, laser and cavity, mutually cancel [3l|, HH, 

m. 
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FIG. 3: Excitation spectrum as a function of the laser detun- 
ing A in the reference frame of the atom. Here, g = O.57, 
k = O.OI7. In (a) A c = -IO7 and (b) A c = 1.2-y. The 
vertical bars indicate the frequency Ao of the carrier (central 
line) and of the red and blue sideband transitions, Ao + v and 
Ao — v, respectively, when the laser is set at A = Ao and the 
trap frequency v — O.27. In (a) Ao ~ A c — v\ In (b) Ao = A c , 
namely cavity mode and laser are resonant. See text. 

We now consider the center-of-mass motion of an atom 
in a harmonic oscillator, and first assume that the me- 
chanical effects are only due to the laser, while the cavity 
wave vector is orthogonal to the motional axis. In this 
regime, the motion gives rise to a modulation of the laser 
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frequency at the trap frequency v. In the regime of strong 
confinement (Lamb-Dicke regime) this gives rise to two 
sidebands of the carrier, i.e. the laser frequency. The car- 
rier and sideband positions are indicated by the vertical 
bars in Fig.[3]in the reference frame of the atom. The cen- 
tral bar is the carrier. The bar at the right (left) of the 
carrier corresponds to a transition which lowers (rises) 
the atomic vibrational excitation by one phonon, namely 
the so-called red (blue) sideband transition. These two 
components are out of phase with respect to the carrier. 
In the limit in which the atomic motion weakly perturbs 
the internal and cavity dynamics, the scattering along 
the sidebands is proportional to the corresponding value 
of the excitation spectrum. Cooling is thus obtained by 
realizing a large gradient between scattering rates along 
the sidebands. Figures [3]Ja) and (b) show two possible 
scenarios, which are discussed in this paper. Case (a) 
corresponds to use the narrow resonance for implement- 
ing sideband cooling with the dressed states [l5|. This 
scenario is obtained by choosing a large value of | A c and 
setting the detuning between the cavity and the laser 
equal to the trap frequency, such that the red sideband 
absorption falls at the center of the narrow resonance. 
This case has been studied in [ll|. In case (b) a large 
gradient is achieved by exploiting the interference profile 
arising when laser and cavity are resonant. 

The dressed state picture, as obtained by neglecting 
the motion of the atom, can be also applied to get some 
insight into the cooling dynamics when cooling is due 
only to the resonator forces or to both laser and res- 
onator. Nevertheless, it does not explain other cooling 
dynamics, which we discuss in this article, and which are 
due to correlations in the gradients of the fields over the 
atomic wave packet. At this purpose one has to consider 
also the quantum motion. 

Figure [4] summarizes the basic scattering processes de- 
termining the cooling dynamics in the basis \g,O c ;n), 
±;n), where n is the number of excitations of the 
center-of-mass harmonic oscillator. The process shown 
in Fig. |3Ua) describes absorption of a laser photon and 
spontaneous emission, whereby the change in the center- 
of-mass state is due to the recoil induced by the sponta- 
neously emitted photon. The scattering rate is scaled by 
the geometric factor a and is found after averaging over 
the solid angle of photon emission into free space. This 
contribution is diffusive, as the motion can be scattered 
into a higher or lower vibrational state with probabilities 
depending on the overlap integrals between the initial 
and the final motional states after a photon recoil. In 
section [V] we discuss the parameter regime in which this 
contribution can be suppressed by an interference effect 
in the dressed states absorption. 

The processes depicted in Figs. 0|b) and (c) describe 
scattering of a laser photon by spontaneous emission 
where the motion is changed by mechanical coupling to 
the laser, (b), and to the cavity(c), field. Since the final 
state of the two scattering processes is the same, they 
interfere. In addition, each term is composed by multi- 



ple excitation paths, and can vanish in some parameter 
regimes. In section [V] we discuss interference effects in 
these two terms. 

The processes depicted in Figs. U^d) an d (e) describe 
scattering of a laser photon by cavity decay, where 
the motion is changed by mechanical coupling to the 
laser, (d), and to the cavity(e), field. The scattered pho- 
ton is transmitted through the cavity mirrors into the 
external modes, and therefore these two processes do not 
interfere with the ones above discussed, but add up co- 
herently with one another. In section IIV CI we discuss 
parameter regimes where interference effects in these two 
terms relevantly affect the dynamics. The general dy- 
namics are a competition of all these processes, and will 
be discussed in detail in the following sections. 

III. THE MODEL 

A. Basic Equations 

We consider an atom of mass M, which is confined by 
a harmonic potential of frequency v inside an optical res- 
onator. The relevant center-of -mass dynamics are along 
the £-axis, while the degrees of freedom of the transverse 
motion have been traced out, assuming that the trans- 
verse confinement is much steeper. Later on we discuss 
how the treatment can be generalized to three dimen- 
sional motion. The atom internal degrees of freedom, 
which are relevant to the dynamics, are the ground state 
\g) and the excited state e), constituting a dipole transi- 
tion at frequency loq and linewidth 7. The dipole couples 
with a cavity mode at frequency lo c and with a laser at 
frequency ujl, whose wave vectors k c and ki, form the 
angle 9 C and 9^, respectively, with the rr-axis. The sys- 
tem is sketched in Fig. [I] and [2] We denote by p the 
density matrix for the atom and the resonator degrees of 
freedom in the reference frame rotating at the laser fre- 
quency. The density matrix p obeys the master equation 

SL p =±[H,p]+C s p + lCp = £.p (4) 

where C is the Liouvillian describing the total dynamics. 
Here, the Hamiltonian H is 

H = H mec + H at + H cav + iJ a t_ ca v + Hl (5) 

where the terms describing the coherent dynamics in ab- 
sence of coupling with the e.m.-field are 

#mcc = |^ + 1 -MV 2 X 2 (6) 

H at = -hA^a (7) 
-H^av = -hSc^a (8) 

Here, x, p are position and momentum of the center 
of mass; a = \g)(e\, and its adjoint; a, a) are the 
annihilation and creation operators of a cavity photon; 
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FIG. 4: Scattering processes leading to a change of the vibrational number by one phonon. The states \g,O c ;n), |±;n) are 
the cavity-atom dressed states at phonon number n. Processes (a),(b),(c) describe scattering of a laser photon by spontaneous 
emission. They prevail in good resonators, for re <C g, 7. Processes (d) and (e) describe scattering of a laser photon by cavity 
decay. They prevail in bad resonators, for 7 <C g, re. The parameters a, tp c , tfiL emerge from the mechanical effects of light and 
arc defined in Sec. UTTbI 



A = u>l — wo an d S c = u>l — u> c are the detunings of 
the laser from the dipole and from the cavity frequency, 
respectively, such that 

A c = A - S c . 

The terms 

-ffat-cav = Tig cos{kx cosQc + <j>) (a) a + aa'*) (9) 
H L = hQ(e ikxcoseL a^ + H.c.) (10) 

describe the radiative couplings of the dipole with the 
cavity mode and the laser, respectively, where g is the 
cavity-mode vacuum Rabi frequency and f2 the Rabi fre- 
quency for the coupling with the laser, is a phase, and k 
is the modulus of the wave vector (|ki,| ~ |k c | m luq/c — 
k). 

The superoperators /C and C s in Eq. ([4]) describe the 
cavity decay and dipole spontaneous emission into the 
modes external to the resonator, respectively, and are 

ICp = — (2apa) — a} ap — pa^ a) (11) 

C sP = 1 (2apo* - o^ap - pa^a) (12) 

where k is the cavity decay rate due to the finite trans- 
mission at the mirrors and 

p= j dcos6 Af(6o)e- ikcose ° x pe ikcoseoX (13) 

describes the events in which the atomic motion recoils 
by emission of a photon at the angle 9q with the trap axis 
with probability JV(8 )d cos 9q. Note that JV(8o) must be 
evaluated taking into account the geometry of the setup. 

For later convenience we introduce the annihilation 
and creation operators b and 6^ of a quantum of vibra- 
tional energy, such that 

x = ^h/2Mv{tf + b) (14) 
p = i^JhMv/2{tf - b), (15) 



and the Hamiltonian term ^ can be rewritten as 

H mcc - hv (tfb + (16) 

We denote by \n) the eigenstates of H mcc at the eigen- 
value (n+ l/2)hv and introduce the Lamb-Dicke param- 
eter 

n = k\-£— (17) 

which scales the mechanical coupling of radiation with 
the atomic motion. 



B. Reduced equation for the center-of-mass 
dynamics in the Lamb-Dicke limit 

We assume the Lamb-Dicke regime, namely the atom 
is localized on a length scale which is much smaller than 
the light wave length, and identify in the Lamb-Dicke 
parameter 77 the perturbative parameter, that allows us 
to treat the coupling of the external degrees of freedom 
with the cavity and the atom internal degrees of freedom 
in perturbation theory [34], [35[ . We apply the formalism 
first applied in [35[ and then further developed in [H, H3, 
l38j ] . Below we summarize some steps. 

At zero order in the Lamb-Dicke parameter the center- 
of-mass is decoupled from the internal and cavity degrees 
of freedom. In fact, denoting by £q = £|,j = o the Liouvil- 
lian at zero order in the expansion, this can be decom- 
posed into the sum of a term acting over the external and 
over the cavity and dipole degrees of freedom, namely 

-Co = £qe + C-oi, 

where 

£oeP = ^r[H mccl p] (18) 
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£■01 P — T^T [H a t + -ffcav + #Oat-cav + Hql, p] 

+K,p + £ 0sP 
and where the Hamiltonian interaction 

-ffoat-cav = ftg(a^<7 + ad 1 ) 

H 0L = /Lfi( CT t + a) 



(19) 



(20) 
(21) 



and the Liouvillian for the atomic spontaneous emission 



£osP = - (^crpa^ - cr^ap - per 1 a) 



(22) 



appear at zero order in the expansion in rj. The term 

g = g cos cj) (23) 

is the zero order atom-cavity coupling strength. 

The spectrum of £q is A = Xj + Xe , where Xj are the 
eigenvalues of £qj and Xe are the eigenvalues of £qe- 
The stationary state is a right eigenstate at eigenvalue 
zero, as it fulfills the secular equation £op — 39]. The 
corresponding cigenspace is spanned, for instance, by the 
eigenvectors p n = pst®\n)(n\, where pst fullfills the equa- 
tion 

£oiPst = 

while the operator |n)(n| is eigenvector of the superop- 
erator £qe at eigenvalue Xe — 0. The corresponding 
eigenspace is infinitely degenerate. We denote by P the 
projection operator over the A = eigenspace, defined as 

Pp = Pst®J2 l">HTr, {(n\p\n)} (24) 

n=0 

where Tij is the trace over the dipole and cavity degrees 
of freedom. At second order in 77 one gets a closed equa- 
tion for the center of mass dynamics of the form 

j t P = V 2 [(S(u) + D)(b^ -fctty) 

+(S(-f) + D)(tfpb - btffi) + H.c] (25) 

where /i = Trj {Pp} is the density matrix for the center- 
of-mass variables, obtained by tracing over the dipole and 
cavity degrees of freedom, and where the coefficients are 
given by 

D = alTni^apst} (26) 
1 I" 00 

S (y) = / dTe iVT Tn {V x e c ° IT V lPss } 
n Jo 

= -TniKiCar + iuyWipss} (27) 
In Eq. ([26]) we used 



a 



d cos #0 cos 8qJV(cos8o), 



which gives the angular dispersion of the atom momen- 
tum due to the spontaneous emission of photons. The 
operator V\ in Eq. ([27]) is given by 



where 



Vi = <plV l + (p c V c 



V L = ifift^t - a) 
V c = —%g{aa* + a/ a) 



(28) 



(29) 
(30) 



describe respectively the mechanical effects of the drive 
and of the cavity at first order in 77, with the two coeffi- 
cients 



tp L = cosV L 

tp c = cos C tan < 



(31) 
(32) 



which depend on the geometry of the setup. Opera- 
tor (|2"8"j) is the gradient of the atom-field interaction at 
the center of the trap and corresponds to the mechanical 
force in the semiclassical limit 401. 



C. Rate Equation 

From Eq. ([25]) one can directly derive the rate equa- 
tion for the occupation probability p n = (n\fi\n) of the 
phonon number state \n), namely 



dt 



where 



Pn = rj 2 [(n + l)A_p. 



Yi+l 



-{{n+l)A + + nA-)p n + nA+p n _ 1 ] (33) 



A± = 2Re{S{Tv) + D} 



(34) 



are the rate of heating (A + ) and cooling The so- 

lution of this type of equation is well known [34]. The 
average phonon number obeys the equation 

(n) = - V 2 (A_-A + )(n)+ V 2 A + (35) 
which, for A- > A + , has solution 

(n) t = {n) e- wt + (n) st (l-e- wt ) (36) 
Here, (n) is the initial average phonon number and 

<->. = m 

is the average phonon number at steady state, while 



W = r] 2 (A- - A 



+ ) 



(38) 



is the cooling rate. 
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D. Discussion 

In Eq. (|33p the internal dynamics enter through the 
coefficients S{v) and D, which determine the rates 
The function S(u) is the spectrum of the fluctuations 
of the radiative force on the atom, namely the Fourier 
transform of the autocorrelation function of the opera- 
tor V\ in Eq. ([28]). For the atom coupled to an optical 
resonator and driven transversally by a laser, we use the 
definition ([28]) in Eq. (|27|). and obtain 



Here. 



St 



-Tn{V L (C 



01 



V L pst} 



is the contribution of the mechanical effect due to the 
laser, the term 



S c = —Tii{V c (C 



01 



i^) 1 V c pst} 



the contribution of the mechanical effect due to the res- 
onator, and 

S c l 



-Tti{V c (Coi + w) 1 V LPSt } 
-TniVLiCoi + i^Vcpst} 



the contribution due to correlations between the mechan- 
ical effects of laser and resonator. Depending on the ge- 
ometry of the setup, one term can be dominant over the 
others. 

The coefficient D, Eq. (|26|) . gives the diffusion in the 
dynamics of the center-of-mass motion. It is the prod- 
uct of two terms: the spontaneous emission rate of the 
excited state into the modes of the e.m. -field, and the 
stationary excited state population, which is determined 
by the overall dynamics at zero order in the Lamb-Dicke 
expansion. 



IV. CAVITY COOLING OF TRAPPED ATOMS 



A. An explicit form of the rate equation for cooling 



An analytical form for the rates entering equation 
can be derived in the limit of a weak laser drive. The 
main steps of the derivation are reported in Appendix A. 
In this limit the heating and cooling rates take the form 

A± = 7a|T s | 2 + 7 |^ L T^ ± + ^> ± | 2 (39) 



with 



0- 



i«;/2 



/(O) 



12* = in^ + Wa) 



(40) 
(41) 



2 L 



7-7, ± 

c 



= in 



= -ci 



g 2 (25 c =F v + i/t) 



g[(AT^ + h/2)(S c + i K /2) 



and 



fix) = (x + S c + vt/2)(x + A + i 7 /2) 



(42) 
(43) 
(44) 

(45) 



The analytic form of Eqs. (|39[) - (|44|) allows one for a more 
transparent reading of these complex dynamics, which 
can be mapped back to the processes shown in Fig. |4] 
The rates are the incoherent sum of three contributions: 
The first term, ja\1s\ 2 , describes a change in the mo- 
tional state by mechanical coupling to the modes exter- 
nal to the cavity, namely by the recoil associated with 
the spontaneous emission of a photon. It corresponds 
to the process depicted in Fig. Ufa) and determines the 
diffusion coefficient through the relation 



D = ja\1 s \ 2 /2. 



(46) 



The second term, ^(px/I^' +Lp c 1~ < ' ± \ 2 , describes scat- 
tering of a laser photon into the external modes by me- 
chanical coupling to the laser (1l) and to the cavity 
(1 C ) field. The two transition amplitudes correspond to 
the processes depicted in Fig. QJb) and (c), respectively. 
They add up coherently and may interfere. Note that 
these processes, together with the diffusive process, are 
dominant for k <C 7. 

The third term, k\(Pl1£' + (p c 1* ,:Sz \ 2 , describes scat- 
tering of a laser photon into the external modes of the 
electromagnetic field by cavity decay. The two ampli- 
tudes, appearing in this term, correspond to the pro- 
cesses depicted in Fig.@Jd) and (e), respectively. Also in 
this case they add up coherently and may interfere. This 
term is dominant for k 7. 

Equation (|39|) contains the basic features of the dy- 
namics of cavity cooling of trapped atoms. It has been de- 
rived (i) in the Lamb-Dicke regime, (ii) assuming that the 
electronic states are bound by the same center-of-mass 
potential, (iii) in the limit in which the laser is a weak 
perturbation to atom and cavity dynamics. Moreover, it 
has been derived for one-dimensional motion. However, 
since at second order in 77 the rate equations for the three 
directions of oscillation decouple in an anisotropic trap, 
it can be generalized to three-dimensional motion as it 
holds for any geometry of the setup. Below we show that 
this equation reproduces a nd g eneralizes results found 
in some particular regimes [l3l Il5|. Moreover. Eq. ([39]) 
allows one for identifying new parameter regimes charac- 
terized by novel dynamics that lead to efficient cooling. 
Some of these dynamics will be presented in section [Vj 



7 



B. Cooling in the bad cavity limit 

Cooling in the bad cavity limit, as discussed in [l3[, is 
recovered by maximizing the ratio A_/A + in the limit 
in which spontaneous emission is negligible. In Eq. (155)) 
we set 7 = and take \5 C \ ^> v. The equivalence with 
the cooling and heating rates reported in [l3[ is evident 
by using the definitions 5 — g 2 A c /(n 2 /A + A 2 ) and 7 = 
<? 2 k/(k 2 /4 + A 2 ), and imposing ip c = lpl = 1 (namely, 
the cavity axes and the laser arc parallel to the atomic 
motion). Below we use this notation but keep tfL and 
ipc thereby allowing for a more general geometry. From 
Eq. 4HD we find 



A± 



a± 



n 2 j 



[(A - <5) 2 + 7 2 /4] ' 

2(6 - i 7 /2) 



(47) 



1 



A-5=F^ + i7/2 
A-6 + i 7 /2 
A-5=F^ + i7/2 



Here, two processes interfere, namely the process in 
which the vibrational number changes by one phonon 
by absorbing a laser photon, depicted in Fig. Hfd) and 
described by the term T£ in Eq. (|39|) . and the process 
in which the vibrational number changes by one phonon 
by scattering a cavity photon, depicted in Fig. (He) and 
described by the term in Eq. ([39]) . These dynamics 
are due to correlations between the mechanical effects of 
laser and the cavity, and they depends critically on the 
geometric setup, as it is visible from Eq. (|4T|) . 



C. Sideband cooling in the good cavity limit 

We consider the case where the atom is far-off reso- 
nance from the cavity and the laser, |A| ^> 7, <?,«;. At 
leading order in A the rates of heating and cooling take 
the form 



" 2 



0(1/A 3 



(48) 



In this parameter regime there is no relevant contribu- 
tion to the mechanical effects from correlations between 
cavity and laser dipole force. Cooling is found for 5 C < 0, 
and the corresponding average phonon number at steady 
state is 



(n) 



(A) _ X + (je + V? 



St 



4(-5> 



with 



D 



g 2 n 



a + (p z L 



+ {5c - vf 



(49) 



(50) 



In the following we do not discuss the solutions leading to 
Dopplcr cooling, and focus onto the parameter regimes 
that lead to ground state cooling, assuming 7 > v. 
For k v Eq. (|49p reaches the minimum value at 

5c = -v, 



where 



I61/ 2 



1 



51) 



(52) 



is the one-atom cooperativity [4l|. The corresponding 
cooling rate is 



W 



|(53) 



Therefore, large ground-state populations and large cool- 
ing rates can be achieved for v 3> k and C\ 3> 1, namely 
for good cavities and in the limit in which the cavity 
linewidth is much smaller than the trap frequency. 

Insight into these results can be found by using the 
dressed state picture discussed in Sec. [TTJ For n <C v, 
large A and 5 C = — v the excitation spectrum corre- 
sponds to the situation depicted in Fig [3)(a) , where the 
the red sideband transition is resonant with the narrow 
resonance at frequency A + = v 1 while the carrier and the 
blue sideband are driven far-off resonance. Hence, this 
condition is analogous to sideband cooling, whereby now 
the narrow resonance is the dressed state of the system 
composed by cavity and atom. 

The cavity loss rate sets the lower limit to the width of 
the narrow resonance, on which sideband cooling is made, 
and thus to the efficiency of the process, as it is visible in 
Eqs. (|5Tj) and (|53|) . From these equations it is also visible 
that large cooperativities ensure better efficiencies. 

The results reported in this section have been obtained 
from rate equation l|3"3"|) with the coefficients Eqs. (|3"9"1) - 
(|44|) expanded at leading order in 1/A. Such expansion 
is valid in the limit where the detuning between atom 
and cavity (respectively, laser) is the largest physical pa- 
rameter. We remark that Eq. (|5ip is in agreement with 
the result reported in [l5| in the corresponding parameter 
regime, whereby the different numerical factors, as well 
as the dependence on the angles, are due to the different 
laser configuration there considered. 



V. COOLING BY QUANTUM CORRELATIONS 
IN THE GOOD CAVITY LIMIT 



In this section we present and discuss novel cooling 
dynamics based on quantum interference effects, which 
are dominant in the good cavity limit, when 

k < 1/, 7, 5. 



8 



In this regime we focus onto the first two terms of the 
rates ([31))) , corresponding to the processes in Fig. 0Ja)- 
(c), and treat cavity decay, giving rise to the processes 
depicted in Fig. 0Jd),(e), as small perturbations. In the 
following we assume that 

<^~(2n + l) 

namely the coupling to the cavity mode does not vanish 
at zero order in the Lamb-Dicke expansion, and it is given 
by Eq. (p?3")) . Comparisons among the efficiencies of the 
cooling schemes are presented in Sec. I VII 



A. Discussion 

Efficient cooling is achieved by maximizing the rate A- 
together with the ratio A-/A+. In this way, ideally one 
maximizes the cooling rate, Eq. (|38|) . and minimizes the 
average number of phonons at steady state, Eq. (|37|) . In 
general, by inspection of Eqs. ([3T?)) - ([4"4")) one can identify a 
strategy for maximizing A-/A+, consisting in identifying 
the parameters such that the heating and/or diffusion 
processes vanish. In this regime, the optimal parameters 
that maximize A^ are found whenever 



Re{/(^)} = 0, 



(54) 



thereby minimizing the denominator of A-. Physically, 
this corresponds to set the red sideband transition at a 
resonance of the atom-cavity system. This strategy is 
effective when the linewidth of the corresponding reso- 
nance is smaller than the trap frequency. Equation (|54[) 
leads to a condition that relates the cavity detuning S c 
with the atom detuning A, namely 



A opt ((5 c ) 



~9 2 



7 k/4 



5 C + v 



(55) 



where we assume fixed couplings and decay rates, hence 
also a fixed cooperativity. For instance, in the case of 
sideband cooling discussed in Sec. lIVCl the optimal cool- 
ing conditions are reached for b c — —v, corresponding to 
the solution of Eq. ([54]) for A — > oo. In this limit, the 
linewidth of the dressed state resonance which is used for 
cooling is infinitely small, and the steady state occupa- 
tion vanishes accordingly. 

Below we discuss various regimes where ground-state 
cooling is efficient and which may be identified, for the 
corresponding parameter regimes, with an approximate 
solution of Eq. (IBT 



Suppression of diffusion by quantum 
interference 



In this section we discuss a cooling scheme based on the 
suppression of diffusion by quantum interference. The 



cooling dynamics are based on the suppression of the car- 
rier transition, and can be understood with the dressed 
state picture. As discussed in Sec. [Til the carrier transi- 
tion can vanish in the regime in which laser and cavity 
arc resonant. The sidebands due to the harmonic mo- 
tion, however, are weak perturbations in opposition of 
phase with respect to the carrier. Thus, they give rise 
to photon scattering with probability given by the cor- 
responding value of the excitation spectrum, depicted in 
Fig. [2Kb) for some parameter regime. The cooling strat- 
egy is thus to enhance the red sideband over the blue 
sideband absorption, thereby profiting of the suppression 
of carrier excitations, and thus of diffusion. This idea is 
reminiscent of cooling mechanisms based on quantum in- 
terference between atomic levels [30|, Wft. whereby here 
the suppression of the carrier transition is due to the 
destructive interference between the laser and the light 
elastically scattered into the resonator by the atom. 

Diffusion is suppressed when 7g = in Eqs. (|31)j) . 
leading to the vanishing of the diffusion coefficient D, 
Eq. (pZFJ)) . From Eq. ([55)) this occurs when 5 C = and, 
ideally, for k = 0. Let us first consider the ideal condition 
of a lossless resonator. In this case, for S c = the steady 
state average phonon number is given by 



[v(y + A) 



3 2 ] 2 + 7 2 // 2 /4 



4^ A (g 



(56) 



Cooling is achieved when either the relations A > and 
g > v, or the relations A < and g < v, are fulfilled. The 
minimum for (n)^ is obtained when A = A opt (<5 c = 0), 
see Eq. (|55|) . For these values the minimum number of 
phonon at steady state is 



7 



A opt (0) 



16A opt (0) 5 

7 V 



16(.g 2 - v 2 ) 2 
with the corresponding cooling rate 



(57) 



W 



A„ pt( or Wi + i 1 ~ l + (4Ao 1 pt (0)/ 7 ) 2 



(58) 



Therefore, ground state cooling, namely 



(0) 



( n )st 



< 1, is 



found when A = A opt and |A| 3> 7, or equivalently for 
\g 2 jv — v\ 3> 7. This condition can be fulfilled (i) when 
v 3> 7 and (ii) when g 2 3> Jf, which is most interesting 
as it can give ground state cooling even when 1/ « 7. 
Below we discuss these two cases in detail. 

Case (i) corresponds to the so-called strong confine- 
ment regime 28j, namely when the linewidth of the 
dipole transition is smaller than the trap frequency. In 
this case sideband cooling is efficient in free space (i.e., 
in absence of the resonator). Like for sideband cooling 
in free space, the implementation inside a cavity res- 
onator requires A = —v, leading to the final occupation 



v(Q) 

/St 



7 2 /I6i/ 2 
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Case (ii) can be fulfilled in the so-called weak con- 
finement regime [28j , namely when the linewidth of the 
dipole transition is much larger than the trap frequency. 
Result (157)) shows that ground state cooling can be effi- 
ciently achieved. This is to our knowledge a novel regime. 
Here, A ~ g 2 /v, so that we can rewrite the cooling limit 
as 

(")^| Aopt(o) «7^ 2 /165 4 - (59) 

We now discuss how the efficiency is modified for 5 C = 
but k finite. At first order in k the diffusion coefficient 
D = 0. In fact, D = 0(k 2 ), being the stationary pop- 
ulation of the excited state of second order in k in this 
regime [32, 33] . At first order in this expansion the steady 
state average phonon number is 



/St. 



V(0) 
/St. 



(1 + F) 



where the term 



2 



2 K - 



<PL<fic 



v A--A- 



<Pc 



(60) 



(61) 



is the correction at first order in k, C\ is the one-atom 
cooperativity defined in Eq. l|5"2")) . and 



A± 



z/ 2 7 



(62) 



[v(v =p A) — g 2 } 2 + v 2 ^ 2 /A 

Cavity decay increases the linewidth of the dressed-state 
resonances, and it thus detrimental. Nevertheless, for 
high cooperativities the result we find in Eq. (|60p ap- 
proaches the result of the ideal case, Eq. (|59"|) . In partic- 
ular, for g 3> v its value at A = A opt (0) takes the simple 
form 

(4? 



(n) 



st 



Aopt(O) 



A opt (o) 8Ci 



(63) 



showing that the corrections scale with 1/C\. 

It must be remarked that the equations presented in 
this section have been obtained from Eqs. (f39|) in the 
limit of weak coupling. Nevertheless, they are also valid 
when the dipole is driven by a saturating laser field. In 
that case, at zero order in the Lamb-Dicke expansion the 
atom is in the ground state and the cavity in a coherent 
state with amplitude C = —Q/g, such that the steady 
state at zero order is 1311 



Post 



\9,Pc)(9,Pc 



(64) 



The derivation of the rate equation of cooling, obtained 
by making no assumption regarding the strength of the 
laser intensity, is reported in appendix [BJ The result 
agrees with the results reported in this section, which 
have been evaluated from Eqs. (|39|) under the assumption 
of weak laser fields. This agreement is not a casuality: In 
fact, when the conditions for this interference effect are 
fulfilled, the atom is driven well below saturation even 
for strong laser and cavity fields, since they mutually 
cancel at the atomic position [33, [HI]. Nevertheless, the 
sideband transitions take place since they are out of phase 
with respect to the carrier. 



C. Suppression of heating by quantum interference 

In this section we discuss a cooling scheme based on 
the suppression of heating processes by quantum in- 
terference. This interference phenomenon is found in 
the good cavity limit, and corresponds to the situation 
in which the heating processes depicted in Fig. 0Jb) 
and (c) cancel out. The corresponding parameters 
are identified in Eq. (|3"9"j) by imposing the condition 
ipi,T^ ,+ + ip c T^ + =0. This condition can be fulfilled for 
instance when ip^ — 0, namely when the laser is orthogo- 
nal to the motional axis and therefore exerts no force, and 
7^r>+ — 0, namely when the transitions to the blue side- 
band induced by the mechanical effects of the resonator 
vanish. Below we discuss this particular case. 

We assume ifiL = 0. Condition T^' + = is fulfilled 
when S c — v/2 and, ideally, k — 0. In this limit A + — 
~fa 1 7s | 2 and the average phonon number at steady state 
is 



(0) 9 7 V/16+ [ff 2 -3t/(A + ^)/2] 
/«/ = a ; — ^ 1 



(65) 



It reaches a minimum for A = A opt (^/2), namely 
A pt(^/2) = 2g 2 /Sv — v, that has the form 



9a 7 2 ^ 2 
A^t (v/2) 16</? 2 16g 4 
with the corresponding cooling rate 



(66) 



W 



Aopt(f/2) 



16f] 2 



,n 2 



9c 



7 (I + 3f 2 /45 2 ) 2 + (3 7 jV85 2 ) 2 

(67) 



To our knowledge, this is a novel cooling regime. Insight 
into these dynamics cannot be simply gained by inspec- 
tion of the excitation rate at zero order in the Lamb-Dickc 
parameter. In fact, the disappearance of the blue side- 
band transition is due to a quantum interference effect 
between the paths of mechanical excitation driven by the 
resonator. Comparing this case with cooling by suppres- 
sion of diffusion, see Sec. lVBl one finds that for S c = v/2 
one can reach lower temperatures in a faster time, as it 
is evident from a comparison of Eqs. (|66|) and (|67|) with 
Eqs. ([IS]), dMD and with the results reported in Sec. UVTH 
A finite, but small, value of k leads to the corrected 
average excitation 



(68) 



where 



F 



9 2 1 



9j 2 v 2 /16 + (g 2 - 3i/(A + v)/2) 
7 2 /4+(A + i^) 2 ) 2(A + i/) 



9 2 1 



71/ 



G 



9 7 2 ^ 2 /16+ [5 2 -3KA + ^)/2f 
Ag 2 ^/v 2 



(69) 



(70) 
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are the corrections at first order in k. They lower the 
efficiency of the mechanism. In particular, the optimal 
final occupation number becomes 



and 



/St 



A„ p tW2) 



1 



557> " 



A op t(f/2) 



<*/<P 2 c + 9 
64Ci 



(71) 



while the corrections to the cooling rate scale with 1/Ci . 
Therefore, for large cooperativities this interference effect 
is relevant to the cooling dynamics. We remark, that as 
for S c — vjl the heating transition vanishes, similarly for 
5 C = —vjl the cooling transition cancels out. 



D. Suppression of diffusion and heating by 
quantum interference 

We finally discuss a cooling scheme based on the sup- 
pression of both diffusion and heating transitions by 
quantum interference. Let us first consider suppression of 
the carrier excitation, that leads to a vanishing diffusion 
coefficient. This can be achieved by using a standing- 
wave drive, such that the trap center is at one of its 
nodes. Therefore, at zero order in the Lamb-Dicke ex- 
pansion the atom does not scatter any photon and the 
cavity is thus empty. Photon scattering originates from 
the dynamics due to the finite size of the atomic wave 
packet, and it is thus a process of second order in the 
Lamb-Dicke expansion. In order to investigate these dy- 
namics, we evaluate the heating and cooling rate entering 
the rate equation (|33|) by considering a new Hamiltonian, 
which is given by operator ([5]) with the new laser-dipole 
coupling 



Hl = fi.fi cos (fcx cos #l + 0i)((7 t + a) 



(72) 



The condition for which the trap is at a node of the laser 
standing wave corresponds to choose 

H = tt/2, 

For this value, the interaction with the laser vanishes at 
zero order in the Lamb-Dicke expansion and the steady 
state is the empty cavity field and the atom in the ground 
state, namely p'§ t = \g,0 c ){g,0 c \. Note that no assump- 
tion has been made on the value of the Rabi frequency Q. 
By expanding Eq. (f?2")) at first order in the Lamb-Dicke 
parameter, we obtain in place of the operator (f2"5|) the 
new interaction term 

V L = -hntp L (a^ + a) 

In the rest of this section we will consider <pL = 1. 

Following the lines of the derivation as in section Ull Bl 
with the new definitions, we obtain the equation for the 
external dynamics, Eq. (|25p . where now 

D = 



S(v)=-Tr I {V L (£' 0I +iv) l V L fJ S t} 



(73) 



Here, C' 0I has the same form as Coi in Eq. with 
Hql = 0. Clearly, the disappearance of the diffusion 
term is due to the fact that there is no field at zero order 
in the Lamb-Dicke expansion. 

The term (|73p . giving the mechanical action on the 
atomic motion, originates solely from scattering of laser 
photons. In fact, the mechanical effects of the resonator 
field appear at higher order in the Lamb-Dicke expansion. 
The cooling and heating rates A'± = 2Re{5'( = F l/ )} take 
the form 



A' 



-y\T£\ 2 + K\Ttf\* 



where now 



7,± 



1L 



z 



K,± 



1L 



n 



+ itt/2) 



and 



f(x) = {x + S c + i/t/2)0 + A + i 7 /2) - g 



-2 



(74) 

(75) 

(76) 
(77) 

(78) 



Hence, the transition amplitudes do not relevantly differ 
from the ones in Eqs. (14"0"1) - (|4"4"]) . However, no low satura- 
tion limit is needed in the derivation of these results. 

We study now the cooling dynamics in an exemplary 
limit, namely in the case of a very good resonator. We 
first consider an ideal resonator, namely n = 0. We ob- 
tain 



A' ± L =0 = 



" 2 (^T^7 m 

f(*T") 2 /4+PTf) (ATf)-9 2 f 



Thus, when 5 C = v then the heating transition vanishes. 
Since the diffusion is also zero, then A', = and 

The corresponding cooling rate reaches the maximum 
value for A = A op t(^), namely A op t(^) = {g 2 — 2v 2 )/2v, 
and takes the form 

W = 4t? 2 — . 

7 

This result is exact for k = 0. 

Finite values of k introduce corrections to the heating 
rate, which takes the form 

7 Ci 

Correspondingly, the average number of phonon at 
steady state is 



(n) 



[2i/(A + v) - g 2 } 2 + 7 2 ^ 2 1 



St 



4i/ 2 7 2 



Ci 



(80) 
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the minimum value of the number of excitations at steady 
state is found at A = A opt (^) and is given by 



These cooling dynamics are novel, and correspond to the 
situation in which the excitation pathways of the com- 
bined multilevel atom-cavity system interfere destruc- 
tively, thereby suppressing the blue sideband excitation. 
They are reminiscent of cooling schemes for multilevel 
atoms discussed in [42| , where suppression of the carrier 
and blue sideband transitions is achieved by quantum in- 
terference between atomic excitations. In the case stud- 
ied here, however, the mechanism which leads to suppres- 
sion of the carrier transition is different from the one that 
leads to the suppression of the blue sideband transition, 
and both are due to the composite effect of cavity and 
laser on the atom. Moreover, the parameter regime here 
considered is one of several possible, that can be identi- 
fied by imposing the disappearance of the blue sidebands 
transition. 



VI. RESULTS 

In this section we compare the cooling efficiencies in 
the various regime, as evaluated from the analytical re- 
sults, and check the range of validity of the analytical 
calculations with a quantum Monte-Carlo wavefunction 
method, where the full quantum dynamics of master 
equation (fj| is simulated. We focus onto the good cavity 
limit, in particular onto the parameters k <C v <C 7. 



A. Plot of the analytical results 

The plots in Figure O show the phonon number at 
steady state and the cooling rate for different geometries 
of the setup. In particular, the plots of the first row de- 
pict the situation in which the mechanical effects on the 
atom are due to both laser and cavity field, the plots of 
the second row show the dynamics when the effects are 
solely due to the resonator, and the ones in the third row 
when the effects are solely due to the laser. The contour 
plots show most evidently the parameter regions where 
cooling is effective. Here, the dashed line represents the 
function which determines the parameters minimiz- 
ing the steady state temperature. Clearly, in the neigh- 
bourhood of this line the lowest temperature is achieved 
in all three cases. We note that the parameter regimes 
where cooling occurs may differ depending on whether 
the dipole forces are due to the resonator or to the laser. 

We now discuss the dynamics in detail. Due to the 
wealth of behaviours, we focus onto the parameter re- 
gions where ground-state cooling appears efficient. 

Figures [5£b),(f), and (j) display the value of the 
average phonon number as a function of 5 C and 



Aopt(^c), namely its value along the function (|55|) . Fig- 
ures GHd),(h), and (1) show the corresponding cooling 
rates. Each plot displays two curves, which have been 
evaluated for two different values of the cavity decay rate 
k (solid curve: k = O.Olv, dashed curve: k = O.lv). 
From these curves it is visible that, as the cavity decay 
rate increases, the cooling efficiency decreases, namely 
the temperature gets higher and the cooling rate lower. 
Nevertheless, for the parameter here considered the cool- 
ing dynamics remains efficient. Let us now discuss the 
behaviour as we vary 5 C and keep A = A opt (<5 c ). 

In all cases the function (n) St exhibits a minimum at 
5 C = —v at very large values of A. This is the sideband 
cooling regime, discussed in Sec. IIVCI The cooling rate 
at these points is very small, since it scales as A -2 , as vis- 
ible from Eq. (f53"| 43]. This cooling scheme exploits the 
dressed states of the system at zero order in the Lamb- 
Dicke expansion, see Sec. HH and its efficiency is thus 
relatively independent of whether the cavity or the laser 
forces determine cooling. Here, the cooling efficiency is 
very sensitive to variation of 5 C , as visible from the con- 
tour plots. This sensitiveness is due to the narrowness 
of the linewidth of the dressed-state resonance which is 
used for cooling the motion. 

The curves in Figures [S^b) ,(f) show also a minimum 
at S c = v/2. This is due to cooling by suppression of 
the resonator's mechanical coupling to the blue sideband 
transition, see Sec. IV CI Clearly, this minimum is more 
enhanced in Fig. [5jf ) , where cooling is solely due to the 
mechanical effects of the cavity, and does not appear in 
Fig. [5]j), where the resonator field has no mechanical 
effects on the atom. The corresponding cooling rate is 
relatively large, being the atom driven close to resonance 
in this regime. An interesting characteristic, emerging 
from the contour plots, is that these cooling dynamics 
are relatively robust to fluctuations of the parameters, 
showing that ground-state cooling is efficient even in the 
regime in which suppression of the heating transition is 
partial. 

The heating region at 5 C = —v/2, appearing in the 
case in which the mechanical effects of the cavity solely 
contribute to cooling, Figs. G3e)-(h), is originated from 
the same interference effects that give rise to cooling at 
the value S c = v/2, and that for S c = —v/2 leads to 
suppression of the cooling transition, see Sec. IV CI 

Another regime where cooling is effective is found at 
6 C = when the mechanical effects arc solely due to the 
laser. This parameter regime is characterized by small 
temperatures and large cooling rates, as visible from 
Figs. E)^i)-(l). This is the regime in which the carrier 
transition is suppressed by an interference phenomenon 
at zero order in the motion, see Secs.lHland lVBI Cooling 
efficiency is robust against fluctuations of the parameters, 
and appear to be relatively stable as the value of k is in- 
creased, as compared with the sideband cooling case, see 

Pig. ED). 

In general, we can conclude that in the range of values 
of S c around the interval [0, v/2], and close to atomic 
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FIG. 5: Average phonon number at steady state (n) St and corresponding cooling rate W, in units of u, in the good cavity limit, 
for k < v < 7 and for three possible geometries: in (a)-(d) (first row) both cavity and laser fields contribute to the cooling: 
Here 9l = 9 C = 7r/4. In (e)-(h) (second row) the mechanical effects of the cavity solely determine cooling: Here 9l = 7r/2 and 
9 C = 7r/4. In (i)-(l) (last row) the mechanical effects of the laser solely determine cooling: Here 9l = 7r/4 and 9 C = 7r/2. The 
contour plots show (n) st and W as a. function of 5 C and A in units of v. The gradation of grey follows the scale where darkest 
region corresponds to the smallest values, the lightest region to the largest values. The corresponding values are reported at the 
bottom of the figure. The heating regions are not coded and explicitly indicated by the label H. The dashed curve appearing 
in each contour plot represents the curve A op t(<S c ), Eq. (|55|l . The parameters are r\ — 0.1, 9l = w/4, n = v, g — Iv, 7 = 10i/, 
K = O.Olf. The other plots display (n)st an d W as a function of 5 C and A opt (<5c), for the same parameters as of the contour 
plot and k = O.Oli/ (solid line), n = Q.lv (dashed line). 



resonance, the cooling efficiency is relatively high. We 
remark that the final temperature is limited by the ratio 
k/v. This is understood in the dressed state picture, as 
the final limit to the narrow dressed state resonance is 
set by the cavity decay rate k. 



B. Numerical simulations 

The curves reported in Fig. [5] have been obtained from 
the analytical equations, which have been evaluated as- 
suming that dipole and cavity dynamics reach the steady 
state on a much faster time scale than the center-of-mass 
dynamics. In Fig. [S] we verify these results by compar- 
ing them with a full quantum Monte Carlo wave function 
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simulation of Master Equation (|4]) . We see that in general 
the analytical predictions are in agreement with the nu- 
merical results for a vast range of parameters, which are 
experimentally accessible. The discrepancies are small 
and are due to parameter regimes where the adiabatic 
evolution is not fulfilled. The discrepancies affect mostly 
the cooling rate, which varies by an overall factor, while 
the final average excitations of the center-of-mass oscil- 
lator are in agreement. 
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FIG. 6: Comparison between the analytical equations, 
Eqs. (|33p and (|39p . and a full quantum Monte Carlo simu- 
lation of Eq. Q. The curves show the evolution of the aver- 
age phonon number as the function of time in units of v~ x , 
the dashed lines correspond to the analytical predictions, the 
solid lines to the quantum Monte Carlo simulation. The pa- 
rameters are 77 = 0.1, 8l — C — (j> = 7r/4, O = v, 7 = 10//, 
k = Q.lv and (a) g = 10//, 5 C = —l.lv; (b) g — 10//, S c — 0; 
(c) g = 10//, 5 C = 0.5//; (d) g = 50//, S c = -1.1//; (e) g = 50//, 
5 C = 0; (f) g = 50//, 8 C = 0.5//; 



the master equation of the system composed by atom and 
cavity, and whose validity is supported by numerical sim- 
ulations taking into account the full quantum dynamics. 
Our analytical results are valid in the Lamb-Dicke regime 
and when the center-of-mass potential is independent of 
the internal state, like in (7I. [25l [26| . 

The equations we derive reproduce the results reported 
in p^Lflql in their specific parameter regimes. Moreover, 
they allow us to identify new parameter regimes where 
the dynamics of the atomic center-of-mass, coupled to a 
laser and a cavity field, results from a non-trivial com- 
petition of the laser and of the resonator dipole forces, 
which can mutually interfere. These interference effects 
are at the basis of novel cooling schemes, which we iden- 
tify in this paper and which allow to reach very large 
ground-state occupations. The corresponding dynamics 
are reminescent of cooling schemes exploi ting interference 
in multi-level atomic transitions [3(| l42l . |44| . 

It must be remarked that the interference effects dis- 
cussed in this work base themselves on the discreteness of 
the spectrum of the mechanical excitations, which is the 
same for the dipolar ground and the excited states. Dy- 
namics will be substantially modified when the external 
potential depends on the internal state, as the spectro- 
scopic properties of the atom are changed and with it the 
scattering cross section. Preliminary considerations are 
in [37| for the case of trapped multilevel atoms, while the 
effects of state-dependent mechanical potentials on the 
cavity field have been discussed in [45| . 

In this paper we have focussed onto the case of dipole 
transitions, which in free space do not allow one for ap- 
plying sideband cooling. Here, ground-state cooling can 
be obtained in good resonators, whereby the final cooling 
efficiency is limited by the cavity decay rate. The bad- 
cavity case is contained in the analytical equations we 
have derived. In this limit the effect of the laser seems 
to be predominant and for a general configuration the 
optimal cooling corresponds to standard cavity sideband 
cooling [46(. This differs strikingly with the good cav- 
ity limit, where correlations between atom and resonator 
can lead to very efficient cooling. We remark that these 
dynamics are largely modified when the drive is set on 
the resonator. 

To conclude, this work considers the cooling dynam- 
ics to the potential ground state of a single atom inside 
a resonator. In the future we will investigate how the 
collective dynamics of several atoms, confined in a res- 
onator, influence the cooling efficiency, and how these 
can be used to prepare quantum states of the system in 
a controlled way. 



VII. CONCLUSIONS 

We have presented an extensive study of the cooling 
dynamics of trapped atoms in optical resonators. Our 
study is based on a rate equation, which we derive from 
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APPENDIX A: LIMIT OF WEAK DRIVING 
FIELD 

We consider the limit when the atom is weakly driven, 
namely when the Rabi frequency f2 is much smaller than 
all the other physical parameters that characterize the 
internal dynamics. We study the dynamics of the motion 
in perturbation theory in second order in and 77, and 
neglect the terms of order rj 4 fi 2 , r/ 2 £l 4 and higher. 

At zero order in f2 the laser field is zero. Hence, at 
steady state the cavity is empty and the atom is in the 
ground state. The state of dipole and cavity is described 
by the density matrix 

Pst = M c )(<7,0 c |, 

which is the solution of the equation £g 7 pg t = 0, with 

(Al) 

+ #0at-cav, P] + TCP + C-OsP 



£qip = £oi 

1 

m 



£"2=0 
-"at + -ffcav 



and is obtained from £, Eq. ([4]), at zero-order in 77 and 
Q. Following the general procedure described in sec- 
tion QTTB] we define the projector Pp = \g,0 c ) (g,0 c \ ® 
J2 n \ n )( n \ r ^ T i{( n \p\ n )} mto tri e eigenspace of the super- 
operator 



no fO 

L -0 — J ~0I 



£ 



0E 



(A2) 



namely the Liouvillian at zero order in 77 and f2, at the 
eigenvalue zero. The effect of the perturbation in 77 and 
fi is described by the equation 

Pp = CpPp 

where we have eliminated the coupling with other sub- 
spaces, thereby obtaining 



n,m=0 y 



2^ r 



2,2 



<m4 



"14 



(A3) 
(A4) 



At the lowest non-vanishing order only the term £2,2 is 
relevant, and 



£ v — r) 2 -=j£-2,2 



(A5) 



a 2 
9 

= ip 2 L C L + y\£ c + <p c ip L £ cL + £ d iff (A6) 

The subscripts L, c and cL label the terms describing 
processes in which the mechanical effect on the atoms 
are due respectively to the laser, the cavity, and to the 



cooperative action of laser and cavity. They have the 
form 



1~L 


— —-tLilLq 


1~.1L 




(A7J 


jC c 


= —P£-ic£q 


1 f fO- 


1 r /»o— 1 /» 

J-olJ-q J-ol 






—PCi c Cq 


1 £ql£q 


L,\ c L,q L.QL 


(A8) 


£ C L 


= PC^Cl- 1 




£qL 






+PC\ c Cq 




1 £\L 






+PC\ c Cq 


1 £il£q 


1 C-0L 


(A9) 



where the terms which trivially vanish have been omitted. 
Here, 



£qlP - 



-t[Hol,p] 



(A10) 



describe the laser-atom interaction at zero order in 77 
where H l is defined in Eq. (|2"Tj) . 



CilP = -r)-[{tf +b)V L ,p] 



(All) 
(A12) 



describes respectively the laser-atom and cavity-atom in- 
teraction at first order in 77 where Vl and V c are defined 
in Eqs. ([2"9|) -([30 |l . The diffusion due to carrier excitation 
is given by 



-Cdiff — PC-2s£-Q 1 £ql£q 1 -^0L 



where 



(A13) 
(A14) 



(A15) 



-(b + tffp-p(b + tf) 2 ] ^ 

is the Liouvillian at second order in 77 for the atomic 
spontaneous emission. 

Tracing over the internal degree of freedom we obtain 
an equation for the density matrix p — Trj {Pp} for the 
center-of-mass variables, 

Tn {drnPp} = r] 2 [D{bptf -tfbp + tfpb-bb^^ + R.c] 
Tri{CjPp} = ^[Sj^ibptf -tfbp) 

+S j {-v)(tfpb-btfp) + R.c.} (A16) 

with j = {L, c, cL}. The coefficients D and Sj are defined 

as 



D = 



T a -r f t r 



0—1 p 

01 Z - 



on£-oi £onP°ss} (-A-17) 



and 



S L (u) = —Tti{Vl (£qi + iv) VlPss} 
S c {v) = -Tr j {V c (£° j + tv)- 1 [V c £° J 1 Col 

+£ol(£q I + ivy 1 V c £^J 1 £ QL p s t} 
S cL (v) = Tyj {v L (£° 0I + ivy 1 Vc^J^oLPlt 

+V C (£° 0I + iv)' 1 £ 0L (£° 0I V LP ° st 
+V C (£01 + if) -1 V^UtPlt} (A18) 
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Setting S{v) = <p1S L {v) + <plS c (v) + (p c (p L S cL (v) we find 
an equation of the same form as Eq. (f2"5)) . The real part 
of these terms are 



Re {D} = ia \T s \ 2 /2 



(A19) 



Re{S L ( T v)} = { 1 \T2' ± \ 2 + k\T^ ± \ 2 )/2 (A20) 
Re{5 c ( T ^)} = ( 7 |T^ ± | 2 + K |T c ^ ± | 2 )/2 (A21) 
Re{S cL (Tv)} = (in^V^* + ^ ± T c ^*)/2 

+c.c. (A22) 

where the coefficients Tj are given explicitly in Eqs. (140|) - 
(144]) . Finally, heating and cooling rates are given by 
y4± = 2Ke{S(^v) + D} and their explicit dependence 
on the physical parameters are reported in Eq. 



APPENDIX B: LIMIT OF SMALL CAVITY LOSS 
RATE 

In this appendix we discuss the derivation of the rate 
equation, in the limit of small cavity loss and 5 C — 0, by 
making no assumption over the laser Rabi frequency f2, 
which may saturate the atomic transition. The results 
we obtain in this appendix are valid provided that g ^ 0. 

In order to derive the rate equation for the atomic 
motion we closely follow the general approach described 
in section 1111 Bl where here we expand at second order 
in r\ and at first order in n. The limit of applications 
of the perturbative expansion are found after identifying 
the smallest rate determining the internal dynamics for 
k = 5 C = 0. This is the width of the narrow resonance 
7_, which for sufficiently large values of |A| takes the 
form 




(Bl) 



Therefore, an expansion in k and r\ is possible provided 
that k <C 7- and rjip c g, rjifLft < 7_. In this regime, 
the internal dynamics at zero order are described by the 
Liouvillian 



0/ 



K=0,5 C =0 



(B2) 



in 



+ H 0L ,p] + C 0s p 



and the steady state, which is solution of -Coo/Post = 0, 
is given by Eq. (|64| . The superoperator at zero order in 
77 and k is given by 



Co 



C 



OQI 



-0E 



(B3) 



and the corresponding projector over the 
eigenspace at eigenvalue zero is Pp = \g, (3 c )(g, (i c \ ® 
^ n \n) (n|Tr/{(n|p|n)}. The dynamics of this sub- 
space at the lowest relevant order in 77 and fc/7- are 



Pp = CpPp where 



E 



n n — C 

7™ 



n— 0,m— 

At lowest order it has the form 



where 



77 2 £2,0 — —PC-iCqqCi 
£-2,1 = PClCnnClCnnlC 



with 



+ PC\Cr ) }KC () }C\ 



C 1 = -r,-[^ + b)V 1 ,p] 



(B4) 

(B5) 

(B6) 
(B7) 

(B8) 



Liouvillian at first order in r\ for the coupling be- 
tween atom and electromagnetic field, and V\ given by 
Eqs. (f2"5|) . Note that we have omitted to write the terms 
which trivially vanish. Tracing over the internal degree 
of freedom we obtain the equation 

A = Trj {Cpposs ® /i} 

for the center-of-mass variables density matrix p = 
Tri {Pp}, whereby 

Tr z {£ v Pp} = if[S t {v){bptf -tfbp) 

+S t (-v){tfph - btfp) + H.c] (B9) 

Here the index I = {0, 1} indicate the order of the ex- 
pansion in k. From this equation we can identify the 
coefficients of Eq. (1231) , and thus 



and 



where 



S (y) 
Si(u) 



D = 

-Trj{Fi (£00/ + i^y 1 Vxposs} 
Tr/jVi (£00/ + 



Pass} 
(B10) 



Note that D = o(k 2 ), as in this system the population of 
the atomic excited state grows quadratically with k (32l . 

Heating and cooling rates are found from the relation 
A± = 2Re{S'o( z R) + SiiTv)}, and take the form 



A ± = tfA ± + + 



(Bll) 
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whereby <C 1. In particular, 

A± = [i/(^A)-j 2 ] 2 |i/ 2 f/4 (B12) 

with 

-t + HF Tl ) Wc (B13) 



Result (jBllI) coincides with the one obtained from ex- 
panding Eq. (|39p to the first order in n and with 5 C = 0, 
from which the results in Eqs. (|58|) and (|60|) have been 
obtained. Nevertheless, in deriving rates (|B11[) we have 
made no assumption on the strength of the laser inten- 
sity. 
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